Strong electron-phonon interaction in multiband superconductors. 
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We discuss the effects of anisotropy on superconducting critical temperature and order parameter 
in a strongly coupled regime. The multiband representation is used as a model for anisotropy. 
We show that strong coupling effects in multiband superconductors lead to pair-breaking due to 
interband coupling because soft phonon modes play the same role as usual impurities. This effect 
makes the order parameters in different bands equal to each other and limits the upper bound on 
critical temperature. 
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INTRODUCTION 

Effects of anisotropy on superconducting critical tem- 
perature and energy gap become of primary importance 
by approaching the strong-coupling regime when transi- 
tion temperature Tc becomes of the order or even larger 
than the characteristic energy 57 of a boson modes which 
mediate superconductivity. This issue received little at- 
tention up to now. In the weak-coupling limit, the effects 
of anisotropy were investigated shortly after the Bardeen, 
Cooper, Schrieffer (BCS) theory ( see, e.g., Ref. [1], and 
for multiband systems Refs. [2], [3]). Following the paper 
by Markovitz and Kadanoff [4], different authors (ref- 
erences can be found in the review [6]) introduced the 
so-called separable interaction 



^^kk' = (l + ak)V'(l + ak'), 



(1) 



where at is an anisotropy parameter, the Fermi surface 
averaging (ak) being equal to zero. The result is the 
enhancement of the effective coupling constant 



Ae// = (iV(0)1/kk'> 



N{0)V{l + {al)) > N{Q)V 



and corresponding rising of the of Tc according to the 
standard BCS expression 



T, = lAAdo exp(-l/Ae//) 



(2) 



where do is the phonon cut-off. 

For multiband clean systems in the weak-coupling limit 
the effective coupling constant in Eq.2 is determined by 
the maximum eigenvalue of the matrix Aq/j, where a,/? 
are band indices. Intraband impurity scattering does not 
affect superconducting properties (Anderson's theorem), 
while the interband one averages out the order param- 
eters Aq and Xeff (and Tc) corresponds to the average 
value 



(A) 



Sa/3^o(0)Aa/3 



(see e.g. Refs. [7], [8]). For positively defined matrix A^^ 
the maximum eigenvalue is bigger than (A) , and we have 
the enhancement of Tc for multiband systems in compar- 
ison with the averaged value independently on the sign 
of the nondiagonal matrix elements which determine the 
anisotropic contribution [9]. 

Recent theoretical studies of superconductivity in the 
two-band superconductor MgB2 [12] and calculations of 
covalent metals as the hypothetical hexagonal LiB and 
boron-doped diamond renewed the interest to the prob- 
lem of an upper bound on superconducting critical tem- 
perature in strongly coupled anisotropic systems. Some 
estimates provide values of A in anisotropic superconduc- 
tors as large as 4 (Ref. [10]) or even 25 (Ref. [11]). 

Let us first remind the result for the strong coupling 
approach to isotropic systems. For the case ft << 2ttTc 
(which can occur for large A) real phonons give the pair- 
breaking contributions to the superconducting pairing as 
well as to the quasiparticle renormalization. The largest 
terms corresponding to pair-breaking and quasiparticle 
damping (see Appendix A) cancel each other ([13, 14]) 
and as the result one arrives the following strong coupling 
expression (see Ref. [15]) 



Tc 



cons 



(4) 



(3) 



where in the simplest approximation const — (27r)^^ ~ 
0.15 (numerical calculations give 0.1827). There are 
interpolation expressions connecting strong- and weak- 
couphng limits (see reviews [6, 16, 17]). 

The authors of Ref. [10] have imposed two possible up- 
per bounds on a maximal critical temperature of multi- 
band superconductors: the lower one is determined by 
the averaged coupling constant (3), while the upper one 
is governed by the maximal (positive) eigenvalue of the 
matrix for the first momentum of the Eliashberg func- 
tions a\f^(ijj)Fajj{uj) 

[X^']^^^Mcf,{l) = 2j^ du;u;al^{Lj)Fcfi{u;) (5) 

( for the Einstein spectrum this value is equal to Aq^JI^). 
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The purpose of this work is to analyze sclfconsistently 
the effects of anisotropy on the upper bound on Tg. We 
show that the low frequency phonons play a role similar 
to intraband and interband static impiirities. The latter 
can lead to the suppression of the anisotropy and as a 
result the upper bound on Tc is determined by the av- 
eraged coupling constant. We consider in a more detail 
the applications to the multiband systems. 



GENERAL DESCRIPTION OF MULTIBAND 
SYSTEMS 



The gap functions {ujn) can be calculated within an 
extension of the Eliashberg formalism to two bands 

Aa {oJn) Za [UJn) = TtT ^ ^ , = , 

/3 \u,,n\<u,c y'^m + ^M'^™) 



■a/3 



w„ 



^^4 + A|K0' 



(6) 
(7) 



where 



/•OO 

Aa/3 = 2 / Wal^ij{w)Fai3 (w) dij/[J^ + [uJm - W„)^], 

Za (wn) are the Migdal renormalization functions and 
u>n = 7rr(2n — 1) and the standard Eliashberg func- 
tions decline the STiperconducting properties and thermo- 

dynamical properties are 



k,k',i/ 
X^('^-'^k-k')> 



(8) 



where a,/3 = {1,2,...}, Na{0) is the partial density of 
states per spin at the Fermi energy, 5^'^, is the electron- 
phonon interaction (EPI) matrix element. Defining 
Xa0 = 2/ u)~^ a'^fj{ij)Fai3{uj)duj , wc obtain the partial 
EPI constants. Values (a;„) enter to the expression 
for the superconducting density of states 



7V(a;) = ^7Va(0)Re. 



(9) 

The Eliashberg functions satisfy the following symme- 
try relations 

N^{0)a%{cj)F^fs{u;) = Np{Q)al^{uj)Fp^{uj). (10) 

For T = Tc we have (we also neglect the Coulomb 
pseudopotential) 



Ac [u^n) Za (Wn) = X! X! 

UJm 



Za {(^n) = i -\ X! X! ^a'jSigrUJm, 
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or 



Aa (w„) 



IH X! X! ^a-iSigmon 



/3 Ulm 



XapAp (ijJm) 



Finally 

A„ (n) p{T,) = "')A/3 (n') , (11) 

/3 n'>l 

where for n,n' > 1 the matrix Ba/3{n, n') has a form (we 
have used the symmetry of the gap function A^ (a;„) = 

Aa (-W„)) 



Baffin, n') = 



Xafjjn- ,,') + X„Jji +n' + 1) 
^(2n- l)(2n' - 1) 

Sff{n), (12) 



where 

S0{n) 



2n- 



n-l 



A;3^(0) + 2 ^ Xf3^{m) 



m=l 



, (13) 



and Aq (n) = A^ (n) /y/2n— 1. Here for the sim- 
plest Einstein spectrum with the frequency Q Xafi{rn) = 
Xa^n^/in^ + 4n^T^m^). The value of Tc is determined 
by the equation 



p{Tc) = 1. 



(14) 



STRONG COUPLING 



The simplest way to estimate p{Tc) for superstrong 
coupling is to put in Eq.ll n = n' = 1. In this case we 
have 



p(Te)+VA, 
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/3 



Xaff 



(27rTe)' + 1 



A, 



(15) 



In the isotropic system Xap = XSaff the last two terms in 
the both sides of the equation cancel each other and we 
have a standard expression for superstrong coupling (see 
Refs.[6],[17]) 



_ ^ 



(16) 
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For the nondiagonal matrix Aq,^ this cancellation does 
not occur and the large Aq/j terms play a role of pair- 
breaking (see Appendix A). 

Let us consider, for the sake of simplicity, the two-band 
system. The solution of Eq. (14) has a form 



Tc,2b = — 



2n 



2C 



with the eigenvector 



{Ai,A2} 



A' + VWT4CD 
2A21C ' 



(17) 



where A = A21A11 + An + 2A12A21 + A12A22 + A22, A' = 
Aii(l + A2i) — A22 (1 + A12), B = A21A11-I-A11 + 2A12A21 + 
A12 A22 + A22 , C = 1 -|- A12 -l- A2i,and D — A12A21 — An A22, 
In this case in the order of 0(1/A) (we suppose An ~ 
A12 ~ A21 A22 ^ A > 1) Ai = A2 and 



1 
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- average (0.4667) 

- max eigenvalue (1 .02) 



N,(0)=1/2N^(0) 




ZTT 

where (A) means averaging over both bands 



(18) 



(A) 



(An + A12) iVi(O) + (A22 + A21) ^2(0) 
iVi(0) + 7V2(0) 



and Na{0) are the partial densities of states. For the 
general (non-Einstein ) spectrum we have (A) il^ =^ 

This means that strong coupling leads to washing out 
effects of anisotropy. 

Similar statements were made in Refs. [18] and 
[19] where the authors have considered the momentum 
dependent interaction. In the former paper the sep- 
arable interaction similar to Eq.l [a'^{uj)F{uj)'\^^, = 
a^{Lo)F{uj)g{p)g{p') with (.g(p)) = 1. They got the re- 
sult that the the expression for Tc in the superstrong limit 
reduces to the isotropic one, while the "pairing potential" 
is proportional to (?(p). This contradicts to more general 
statement in the latter article, where the positive (at- 
tractive) interactions a^(k, k', w)i^(k, k', (jj) for all k, k' 
was investigated and it was shown that the gap function 
becomes the k -independent that leads to the isotropic 
expression for Tc- The detail inspection the situation in 
Ref. [18] also shows that the real order parameter which 
enters to the density of states (see Eq.(9)) is isotropic for 
large A. 

We have investigated numerically the evolution of Tc 
and eigenvectors Aq, as functions of the coupling strength 
Aq,^ for the model matrix of the Eliashberg functions 



a^^(w)F„^(w) = a'^{LL!)F{uj) 



1 1/5 
1/10 



(19) 



FIG. 1: (color online) Critical temperature (a) and the gap 
ratio (c) in the two-band case as a function of intraband cou- 
pling constant A in the first band. The panel (b) shows that 
Tc in the strongly coupled regime is determined by the av- 
erage coupling constant. The panel (c) shows the ratio of 
the order parameters in the two bands. The gap function 
becomes isotropic in the strongly coupled regime. The nu- 
merically calculated ratio A1/A2 is very accurately described 
by the expression Eq.(17) in a broad range of A. 



We suppose, for simphcity, 2iVi(0) = N2{0) and A22 = 
(i.e. no intrinsic superconductivity in the second band). 
The average EPI constant (A) is equal to 0.467A. Re- 
sults for Tc are presented in Fig.l. We see that for weak 
and intermediate coupling there is an enhancement of 
Tc due to anisotropic effects in comparison with aver- 
aged value. For small EPI the result coinsidc with the 
weak-coupling expression for Ac// — Amax o.v. = 1.02A, 
where Amax o.v. is a maximal eigenvalue of the matrix 
(19). This enhancement, however, vanishes for large val- 
ues of A when phonons lead to isotropization of the su- 
perconducting order parameter. 

We have to note that the result (18) is obtained un- 
der the condition of nonvanishing (A) and in the Born 
approximation [21] for the spin-independent interaction. 

Recently the model for the system with strong 
coupling anisotropic interaction was considered in 
Ref. [22]. It was supposed that the difference 
between the interaction in the quasiparticle channel 
^Q;^(k, k', w)T'(k, k', Ci;))^^ and in the Cooper channel 
(A(k)a2(k, k', Lu)F{k, k',a;)A(k'))^<j / (A2(k))^^ is in- 
dependent on the coupling strength. Above analyzes (as 
well as [18, 19]) shows that this difference vanishes for 
strong coupling. This removes unphysical results for Tc 
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obtained in this limit in the mentioned paper. 

In Appendix B the sensitivity of to different phonon 
modes is considered by calculating the variational deriva- 
tives. It is shown that the negative (divergent at small 
frequencies) contribution to the nondiagonal variational 
derivative of Tc vanishes in the strongly coupled regime. 

CONCLUSIONS 



Now let us separate the functions a^^(f2)FQ,^(0) on to 

low energy part ^a^^(0)FQ,/3(f2)^ and the high-energy 
one 

The same procedure can be done for the kernel (20) 



We have shown that strong coupling effects in the 
multiband superconductors lead to the appearance of the 
strong damping which results from pair-breaking due to 
interband coupling. 

For systems with the attractive interaction this ef- 
fect leads to averaging of order parameters in different 
bands. As a result asymptotic behavior of Tc is de- 
scribed by the well known single-band expression Tc oc 

= ^2 J^duLO (a'^{uj)F{uj)). This means that 
the upper bound on Tc in the superstrong coupling regime 
is determined by the averaged coupling constant, while 
the higher upper bound corresponding to the maximal 
eigenvalue of the matrix [Afi^] is never reached. 
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APPENDIX A 

We extend the results of Ref. [14] for effects of low 
frequency intermediate boson modes { ^ ^ 2TrTc) on the 
critical temperature of the multi-band superconductors. 
On the real frequency axis the equations for the complex 
order parameter Aa{u!) and the renormalization function 
Za (w) have forms ( we neglect the Coulomb contribution) 

/CO A / /\ 

dzKc.f,iz',uj)Re^^, (Al) 
-oo ^ 



/oo 
dzKaf3{z',L0), 

P "°° 



(A2) 



where Kap{z' .w) is a kernel of the interelectron interac- 
tion via intermediate bosons with the spectral function 



alJn)Fc,p{a) 



1 f°° 

K^0{z',u;) = - / dnalp{n)Fc.p 



(Q) X 



tanh ^ + coth ^ 
z' + — CO — iS 



K„0{z',uj) = K<{z',io) + K>{z',io). 



(A3) 



In the first term in the right hand side of Eq. (A3) we can 
neglect the frequency in the denominator. In this case 



TT Z' 



i5' 



(A4) 



where 



/>OG 

r<^ = TT / dn {al^{n)Fcis{n))^ coth{n/2Tc) ~ 27rA<^Te 

(A5) 

is the matrix of the electron scattering on the low-energy 
excitations. Now we use the dispersion relation for the 
order parameter A^(w) 



dz' 



u! — z' + id 



Re 



(A6) 



which is a consequence of the dispersion relation for the 
electron Green function in the Nambu representation. 
Combining expressions (A1,A2) with (A3-A6) we get 



Aa(w) 



j:'^A,iu;)+Y: r dzK>,[z'M^^ 



We see that the low frequency excitations play a role 
of intraband and interband static impurities. Intraband 
ones drop out from the Eq. (20) (so called Ander- 
son's theorem). It is interesting to note that the famous 
cancellation of the largest terms proportional to A< (see 
e.g., [6]) comes not from the strong renormalization of 
the quasiparticle energy (ReZ), but from the damping 
iV<. 



APPENDIX B 

In Ref. [23] the sensitivity of Tc to different phonon 
modes was considered by calculating the variational 
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derivatives 5Tc/5\a^{^)F{Vl)\^^. For the diagonal el- 
ements (a = (5) the result for small f2 (f2 <C ^ttT^) co- 
incides with the one obtained by Bergmann and Rainer 
[24] 5Tc/5 [a^{n)F{n)\ ~ for the isotropic single band 
system. This corresponds to the enhancement of by 
adding low frequency phonons (bosons). 

In the multiband case the interband derivative has the 
following form 



6T, 



E 

n>l 



+ 0{n), 



and dTj6 [Q^{n)Fin)] and 6TJ6 [a'^{n)F{n)]^^ have 
different signs. This contradicts to the symmetry rela- 
tion (10). If we change the function [a^(ri)F(f2)] the 
counterpart [Q:^(n)F(0)] has to be changed automat- 
ically. Only the symmetrized off-diagonal combination 
dTc/ da^ {n)F{U)o,d. has physical meaning. Here 

a\n)F{n)oA. = 

iVi(O) [a%n)F{Q)]^^+N,{0) [a%n)F{n)]^^ 
7Vi(0) + 7V2(0) 

As a result, we obtain 

ST, ^ 

n>l 



Sa^n)F{n)o.d. 



7V»(0) ^ [Aa (t^n)-Ai {ujr,)Y 



o{n). 



(Bl) 



In contrast to the single band case (see [17]), the off- 
diagonal derivative has different behavior in the weak- 
coupling and strong-coupling regimes. For the former 
case one can suppose A^ (a>„) = Aq(0— |ci;„|) and 

5Tc/5a^{9.)F{VL)oA. -^^^TF^- This means that the 

addition of nondiagonal interaction with low frequency 
phonons leads to strong suppression of the critical tem- 
perature in weak-coupling anisotropic superconductors. 
This result was obtained in Ref. [25] for the anisotropic 
separable interaction. In the strong coupling limit, as it 
was shown above, Ai A2, then the first term in (Bl) 
vanishes and STc/ 5a^{Vl)F{Vt)o.d. ~ O > 0, similar to 
the intraband contribution. This result can be directly 
obtained from the expression (18). 
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